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Abstract 

The linear and the nonlinear non-Markovian quantum state diffusion equation 
(NMQSD) are well known tools for the description of certain non-Markovian open 
quantum systems. In this work, we systematically investigate whether the normalized 
linear NMQSD or the nonlinear NMQSD solutions can be generated by means of a 
time-continuous measurement. By considering any conceivable measurement scheme 
in the framework of instruments, we derive a necessary criterion for a measurement 
interpretation of both equations. Concrete examples show that the normalized linear 
NMQSD solutions are realizable only in the Markovian limit in general. The applica- 
tion of the presented criterion to the nonlinear NMQSD remains an open issue. 
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1 Introduction 



There is no conceptual problem to assign a pure state to a classical open system like a 
Brownian particle. In quantum mechanics however, the entanglement between the open 
system and its environment complicates things. Even if the state of total quantum system 
is pure, the only state which can be assigned to an open quantum system is a mixture: 
the reduced density operator p re d(t). Obviously, the reduced density operator can be 
expressed as a weighted sum over pure states, which is the starting point for powerful 
Monte Carlo methods. As a matter of fact, this so-called unravelling is not unique pQ. On 
the other hand, a measurement of the environmental degrees of freedom can destroy the 
entanglement between the system and its environment. The open system is therefore in a 
well defined pure state directly after the measurement. 

In order to describe the dynamics of a Markovian open quantum system by means of 
Monte Carlo methods, stochastic Schrodinger equations are often employed. The solutions 
of these evolutions equations for pure states have the remarkable property that their 
expectation value coincides with p Te d{t). Yet such stochastic Schrodinger equations do 
not only occur in the context of Markovian open quantum systems but also as wave 
function collapse models and for the description of certain time-continuous measurements 
[2HS]. Moreover, there is a close relation between time-continuous measurements and the 
stochastic Schrodinger equations for the description of Markovian open quantum systems: 
The solutions of the latter evolution equations can be interpreted as trajectories of post- 
measurement states conditioned on a time-continuous measurement. 

In the past fifteen years, stochastic Schrodinger equations have been constructed in 
order to deal with general non-Markovian open quantum systems as well. Yet the under- 
standing of these equations in the context of time-continuous measurements has not been 
established. Before we come to this point in more detail, some preliminaries are presented: 
Firstly, we review a certain stochastic Schrodinger equation approach to non-Markovian 
open quantum systems in section [TTT1 Afterwards, the discussion so far about the relation 
between these equations and measurement theory is presented (section 1 1 . 2[> . In section 
11.31 a general framework for the description of time-continuous measurements is outlined, 
which constitutes the basis for our systematic investigations of the interpretation of this 
particular stochastic Schrodinger equation in sections El and [H 

1.1 Stochastic Schrodinger equations and open quantum systems 

In many settings, some system of interest, characterized by some Hamiltonian H, is cou- 
pled to an environment consisting of finitely or infinitely many harmonic oscillators, where 
the coupling is assumed to be linear in the environmental annihilation and creation oper- 
ators a\ and at, respectively [7]. Here the structure of the environment is determined by 
the oscillator frequencies uj\ and their distribution. The coupling is characterized by the 
coupling strengths g\ which quantify how strongly each environmental mode A couples to 
the open system via a system coupling operator denoted by L. In the interaction picture 
with respect to the environmental Hamiltonian H env = ^2 x w a^1<2a> the total Hamiltonian 
reads 

H tot (t) = H®t + Y,{a*xL® ajy^* + g x L ] ® axe'^ ) (1) 

A 

for this class of open quantum systems. In the following, we consider only scenarios with 
the initial condition for the total system: (^(O)) = \ipo) (£) x \0x)- Here |0a) denotes the 
vacuum state with respect to ax- The environment is hence assumed to be initially at zero 
temperature and uncorrelated with the open system. It turns out that the embedding of 
the system into its environment is characterized by the environmental correlation function 
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a(t — s) which equals 

a(t-s) = Y,\9x\ 2 e- i ^ (t - s) (2) 

A 

for the above initial condition. Note that this correlation function is hermitian, which 
means a*(s — t) = a(t — s). 

If the Born-Markov approximation is applicable, p re d(t) evolves according to a master 
equation of the Lindblad form, where L turns out to be the corresponding Lindbladian 
[3"|I71IB] . In this paper, an open quantum system is called Markovian if and only if its reduced 
density operator obeys a Lindblad master equation with time-independent Lindbladians. 
Such a Markovian open system dynamics can only arise in the microscopic model (HD if 
one formally assumes a(t — s) = n5(t — s). For large system Hilbert space dimensions, 
it is, however, more efficient to employ so-called Monte Carlo wave function techniques 
or quantum trajectories methods instead of solving the known operator-valued Lindblad 
equation numerically [21[3]. These methods are based on stochastic evolution equations 
for pure system states as sketched above. 

The dynamics of an open system becomes non-Markovian if the environment is struc- 
tured and there is no time-scale separation or if the coupling is strong. Measures for the 
non-Markovianity of open system dynamics are currently discussed [9j. In this regime, 
the form of the master equation is in general unknown. Super-operator techniques are 
available for determining p re d(t) but often not practical |10j . Monte Carlo wave function 
methods, however, provide an alternative means for coping with this problem |li y i2 j . In 
this paper, we concentrate on an exact generalization of the Markovian quantum state 
diffusion equations to the general setting (JTJ) which has been found by Diosi, Gisin and 
Strunz [13pi4j. The so-called linear non-Markovian quantum state diffusion equation (lin- 
ear NMQSD) reads for an initially zero temperature environment 

dtMPU = ( - i# + zfL - D J* ds a(t - s)^) \M[z%)), (3) 

where [z*]q indicates a functional dependence on the drawn path of the complex Gaussian 
process z* . This stochastic process is completely characterized by the moments E[zj] = 
E[z(Z s ] = and E[ztZg] = a(t — s). By solving ([3]) for a deterministic initial condition \vpo) 
and many realizations of the stochastic process z*, one can unravel the reduced density 
operator of the open system corresponding to /5 rc( j(0) = \ipo){ipo\- 

^ ed (t)=E[\^%))(M[z*] t o)\]. (4) 

Here E[-] denotes the expectation value over all solutions of ([3]). When reviewing the 
derivation of ([3]), one realizes that for any fixed time t > IV^G^o)) ^ s the unnormalized 
relative state of the system corresponding to a coherent environmental state [13]: Let \z\) 
denote the unnormalized Bargmann state with (z\\z\) = exp(— \z\\ 2 ) corresponding to a\ 
and let d 2 z\ be defined as dD\e(z\)d3m(z\) . The state of the total system at time t can 
then be represented as 

l*M>=/n( d2zAt T ( ~ k|2) ) hMM8»®M P) 

A A 

if [z*]q is given by z* = -iYy\9*\ z \ eluJxT where r € [0,t]. 

In the special case of a Markovian open quantum system, i.e. a(t — s) = n5{t — s), 
the linear NMQSD reduces to the well known linear Markovian quantum state diffusion 
equation driven by white noise [2] . It has to be pointed out that the open system dynamics 
can be non-Markovian according to the definition given previously while z* is a Markovian 
stochastic process. A correlation function for which this is the case will be discussed in 
section 13.21 
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The occurrence of the functional derivative in ([3j) is a direct consequence of the driving 
with coloured Gaussian noise [15]. Yet if there is a system operator 0(t,s, [z*]q) which 
removes the non-locality of ([3]) with respect to the stochastic process, i.e. 

±J M z%)) = d(t, Sl [z%)\M[z%)), (6) 

([SlH]) can be used as a Monte Carlo scheme for calculating p Te d(t). If such a replacement 
is feasible, the linear NMQSD becomes 

dt\A([z*\o)) = ( " + & ~ # f ds a(t - s) 6(t, s, [**]<)) |M*1o)>, (7) 

Jo 

and one can apply parallel computing in order to solve it for many process realizations. 
There are many examples where the replacement ([6]) exists and an exact ansatz operator 
can be stated. Two of them will be reviewed in section [3j In other cases, approximation 
schemes have to be employed in order to calculate 0(t,s, [z*]q) [16J. As we will see in 
section O ([7J) leads to a reversible dynamics of the relative states |V't([- 2 *]o))- This does 
not imply that p ie d(t) evolves reversibly, of course. Since the right hand side of ([7|) is 
not anti-hermitian, the squared norm of a solution iV'tQ^Jo)) is n °t preserved under the 
dynamics. Yet its expectation value remains constant [14J. This observation becomes 
important for the following idea. 

For estimating the necessary sample size for a good Monte Carlo simulation of p re d (t) , it 
is instructive to average over normalized states |V't([- 2; *]o)) = lV'i([ 2; *]o))/ll' ( M[ 2; *]u)ll i n ©• 
If Vq{-) denotes the Gaussian probability measure of the paths [z*]&, flD 

can equivalently 

be written as 

pUt) = /|&M)X&([**]S)I ll^(rt)ll 2 ^(«). (8) 

Here and from now on, we suppress a possible dependence on the complex-conjugated 
process realization [z] in the notation^. According to the observations lined out above, 
the probability weight with which |^t([^*]o)) contributes in i.e. ||V't([ 2: *]o)ll 2l/ o(d[2 ; *]o)) 
changes in the course of time. Consequently, the necessary sample size for a good estimate 
of Pred(t) depends on the instant t, which makes a simulation of the dynamics of p Te d(t) 
difficult. 

The numerical efficiency can significantly be improved by a Girsanov transformation, 
which results in the norm-preserving nonlinear NMQSD [J3]. The latter refers to the 
evolution equation for the normalized states 

l&([*To)> := \M?U, (9) 
where Zt denotes the shifted process 

zf-=z t + f dsa(t - s) (M^r )\m([z*}o))- (10) 
Jo 

This process has to be regarded as a functional of the coloured noise Zt, which we omit in 
the notation. In this sense, the reduced density operator is also unravelled by the nonlinear 
NMQSD solutions: 

Pred(t) = /|^([r]»(rr )| *g(drt). (11) 

Consequently, the Girsanov transformation leads to a stochastic Schrodinger equation 
whose solution occur with a time-independent weight in the unravelling. 

1 The path probability measure Uq actually depends on both [z*]o an d [z]q. 
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In the following, we refer to both the normalized solutions |V>t([<z*]g)) of the linear 
NMQSD and to the nonlinear NMQSD solutions |^t([£*]o)) as (non-Markovian) quantum 
trajectories. As stated above, the expectation values of those objects have a well de- 
fined physical meaning. This paper is devoted to the question whether a single quantum 
trajectory is physically realizable, too. 

1.2 Stochastic Schrodinger equations and time-continuous measurements 

Whether quantum trajectories are subjectively real [IT], i.e. can be interpreted as trajec- 
tories of post-measurement states of a time-continuous measurement scheme, is a question 
motivated by both conceptual and practical interest: The variety of distinct unravellings 
could possibly be understood as a consequence of distinct measurement schemes [151 [TU], 
From an experimental point of view, a measurement scheme which generates quantum tra- 
jectories would feature a non-demolition property: If applied non-selectively upon many 
identical copies of the system, this scheme generates the state of the unmeasured open 
quantum system, i.e. p ie d(t). Certainly, such a measurement scheme would depend on 
a(t—s). So one could simulate the embedding of an open system into differently structured 
environments. 

In the Markovian regime, the interpretation of quantum trajectories is well established: 
Quantum state diffusion driven by real (complex) white noise can be realized by a homo- 
dyne (heterodyne) detection of the environmental photon field [SIS]. There are many 
other examples in which Markovian stochastic Schodinger equations have been postulated 
in order to model imprecise time-continuous measurements or the dynamics of the wave 
function reduction (e.g. [5]). Moreover, Wiseman and Diosi have completely parameterized 
all diffusive stochastic Schrodinger equations for Markovian open systems and linked those 
equations with generating measurement schemes [19] . Some of those equations generate a 
non-Markovian system state dynamics while the expectation value of those states obeys 
a Lindblad master equation. Furthermore, the general form of a stochastic Schrodinger 
equation driven by independent Wiener processes has been derived such that it can be 
interpreted as an evolution equation for post-measurement states [20J . If one averages the 
post-measurement states over all measurement records, the resulting a priori state evolves 
according to a master equation of the Lindblad type. 

However, the most general class of diffusive stochastic Schrodinger equations which 
allow a measurement interpretation has not been characterized in the non-Markovian 
regime. There are examples where such an interpretation is indeed feasible |21| and, 
recently, Yang, Miao and Chen proposed an interesting approach relating non-Markovian 
open quantum systems to time-continuous measurements |22j. In this context, see also [23 . 
Yet in particular, the meaning of the non-Markovian quantum trajectories is a vividly 
debated open question. 

From ([5]), it becomes clear that for any time t > the open system can be forced 
into a pure state proportional to |V ; t([- 2; *]o)) by projectively measuring in which coherent 
state \z\) the environment is. In this context, see also the work of Gambetta and 
Wiseman |24j . who furthermore show that such a single-shot measurement interpretation 
leads to alternative unravellings as well. If, however, those single-shot measurements 
are used to constitute a time-continuous measurement scheme performed upon a single 
copy of the system, the system dynamics will not be given by a NMQSD because of the 
incompatibility of the single-shot measurements. This fact has motivated an interpretation 
of NMQSDs as a hidden variable theory [25J. Furthermore, Diosi suggested that a certain 
readout schedule of initially entangled von Neumann detectors could induce the NMQSD 
dynamics |26j . In response, Gambetta and Wiseman showed that this scheme does not 
preserve the purity of the system states |27j . 

In [24,25,27], Gambetta and Wiseman argue that pure-state quantum trajectories do 
not exist for general non-Markovian systems within the framework of standard quantum 
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mechanics. Their conclusion, however, is based only on a concrete measurement scheme 
which fails to generate non-Markovian quantum trajectories. In particular, the line of 
argument in [24] neglects the possibility of a measurement scheme which adapts condi- 
tioned on previous outcomes. Therefore, the central aim of this work is to investigate 
the realizability of non-Markovian quantum trajectories systematically by considering any 
conceivable time-continuous measurement scheme. 



1.3 General description of time-continuous measurements 

The general framework for describing quantum measurements was studied by Davies and 
Lewis [28], who have introduced the concept of (completely positive) instruments. In 
particular, this concept covers von Neumann measurements |29j but also indirect and 
destructive measurements (e.g. photodetection). A more recent introduction into general 
quantum measurements and instruments can be found in [30J and [31], respectively. 

Let Q be the set of all possible measurement outcomes and let J 7 be a cr-algebra over 
f2 such that it contains all events which can be verified or falsified by a measurement. 
An event F G F is called verified in an actual measurement if the measurement outcome 
uj G is an element of F. Then a map y (•)[•] : J~ x £>(H) — > £>(H), where £>(H) denotes 
all bounded linear operators on the system Hilbert space H, is called an instrument if 
the following conditions are fulfilled: For any F G F, Y(F)[-] is an operation, i.e. linear, 
trace-decreasing and completely positive. Besides, the normalization tr(y(fJ)[A|) = tr(A) 
is satisfied for any A. Finally, y(-)L4] is cr-additive for any A. 

Given the pre-measurement state p, the probability for verifying F G F is postulated 
to be 

P(F\p) := tv(Y(F)[p]) = tr (pY^F)[l]), (12) 

where y^(i ? )[-] denotes the dual map of y(i ? )[-]. The corresponding post-measurement 
state is then 

._ Y(F)[p] 

PF - wmKr (13) 

A dilation theorem allows the interpretation of the stochastic dynamics (|12I13|) as being 
induced by a projective measurement upon a measurement apparatus which has been 
entangled with the system by some interaction [32] . Such an extension is called a measuring 
process and is not in general unique. 

A general time-continuous measurement can be described either by explicitly mod- 
elling the interplay of general single-shot measurements and short periods of free evo- 
lution (e.g. [6]) or as measurements that formally follow the mathematical framework of 
single-shot measurements (e.g. [20134]). Although the latter description of time-continuous 
measurements appears to be rather formal, it is the appropriate means to deal with the 
issue of this paper. The advantage of the latter description in comparison to the former 
is that one need not know which physical quantity has to be measured after the short 
periods of free evolution. 

The chosen approach starts with the continuum limit: The measurement record is 
assumed to be a function of time instead of single numbers from the beginning. So for 
describing a measurement within the time interval [0,t], the set of all measurement out- 
comes needs to be an appropriate space of functions, e.g. all real- or complex-valued 
right continuous functions on [0,t] with left limits. We denote a single record by [x]q G Qq, 
i.e. [x]q := (t i-> x(t),t £ [0, i]). Naturally, the measurement record x(-) must be a non- 
anticipating stochastic process, i.e. adapted with respect to some filtration {Ft}t>o- 

Then a family of instruments y = {Yq(-)[) : F t x £>(H) ->■ £>(H); t > 0} describes the 
whole measurement. In this approach, causality has to be ensured manually: Let Fq G F s 
be some event which can be verified or falsified by a measurement up to time s. The 
probability for verifying Fq must not depend on whether or not one measures longer than 
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necessary, say up to time t > s, and then throws the additional data away. In other words, 
y must fulfil 



P^\p) = tr(Y s (F *)[p]) = P*(F£ x nl\p) = tr(y '(^ x ft*)[fl) (14) 

for any pre-measurement state p. Here Fq x 0* refers to all elements of f2g which lie in Fq 
if restricted to [0, s]. Consequently, y can only represent a time-continuous measurement 
scheme if the following compatibility demand is satisfied 

V0<s<i, VFSeT,: (y s ) t (F s )[l] = (y *) t (F s x^)[l]. (15) 

According to a theorem proven by Loubenets, any instrument can be represented in a 
particular integral form [33J. As the realizability of certain pure post-measurement state 
trajectories is studied in this work, we only state the representation of instruments cor- 
responding to efficient measurements. In this context, a measurement is called efficient if 
a pure pre-measurement state is turned into a pure post-measurement state as a conse- 
quence of the measurement. So if y represents an efficient time-continuous measurement, 
there are operator-valued functions [x]q h-> V^QxJq), the stochastic evolution operators, 
and finite scalar measures pf Q : Tt — > I C [0, oo) for any < s < t such that 

YZ(H)IP]= [ ^(Wo)/5(^(Mo)) f ^(d[x]* ). (16) 
Jf* 

The probability for finding a measurement record in the vicinity of [x]q reads 

P '(d[x] t |^) = ||y *([x]' )^|| 2 /4(d[x]* ), (17) 

given the initial pre-measurement state The corresponding post-measurement state 
equals 

mxU _ . (i 8 ) 



By stating conditions on V^*([x]q) and such that the compatibility condition (fT5|) 

is satisfied, Barndorff-Nielsen and Loubenets have presented a general framework for the 
description of time-continuous measurements [33] ■ When applying this formalism to our 
problem, we start with determining the ingredients of (|16p and then rederive those con- 
ditions on the ingredients from the compatibility demand (|15p . In particular, those con- 
ditions ensure that the instrument for an observation within [0,t] can be represented as 
a composition of arbitrarily many instruments describing a sequence of measurements up 
to time t. Namely for any s € (0,t), one can write 

Y$(d[x]i x d[x]l)[p] = Y s \d[x]l\[xUY s (d[x] s )[p}], (19) 

where ^(■|[ x ]o)N denotes an instrument which depends on the previously measured signal 

Y*(d[x]l\[x] s )[C] := V*([x]l) C /i*(d[x]*|[x]g). (20) 

Here the conditional measure /i*(d[x]*| [x]q) = /J (d[x]Q x d[x]*)//iQ(d[x]o) has been intro- 
duced. The possible dependence of Y* on the result [x]q can be seen as a generalization of 
instrumental processes with independent increments |31|.I35J : The time-continuous mea- 
surements described in [34J are allowed to adapt conditioned on previous measurement 
outcomes. 
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2 A necessary criterion for a measurement interpretation 



In the following, we consider only NMQSDs whose ansatz operator is noise-independent^, 
i.e. 0(t, s, [z*]q) = 0(t, s) with 0(s, s) = L. In order to study the transition from the non- 
Markovian to the Markovian regime, we do not specify the correlation function a(t — s). 

Suppose there is an interpretation of the normalized solution of ([7]) in terms of some 
time-continuous measurement scheme. Then there must be a family of instruments y such 
that the following demands are satisfied: 

i. These instruments are compatible with respect to time, i.e. obey (|15p . 

ii. Given some fixed pure initial pre-measurement state, the set of all post-measurement 
state trajectories coincides with the set of all normalized linear NMQSD solutions 
corresponding to the given initial condition. 

iii. If an event Fq £ Tt is associated with a zero probability set of normalized linear 
NMQSD solutions, Fq shall be verified with zero probability in an experiment. 

The minimality demand (iii) ensures that the measurement scheme generates only those 
normalized NMQSD solutions as post-measurement state trajectories which actually con- 
tribute to the unravelling of the reduced density operator (|8j) . Alternatively we could have 
demanded that the a priori state of the measurement scheme equals the reduced density 
operator. Yet such a stronger demand would blur a result which we anticipate here: Up 
to a certain degree, it turns out to be unimportant whether one investigates a possible 
measurement interpretation of the linear or of the nonlinear NMQSD. Both issues can be 
handled within a single framework - at least in principle. 

The derivation of a necessary criterion for the measurement interpretation of non- 
Markovian quantum trajectories consists of three steps: First of all, it has to be proven 
which quantity can be assumed to be the measurement signal w.l.o.g. (section 12. ip . Sec- 
ondly, the ingredients of the integral representation (|16p are derived (sections 12.21 and !2.3p . 
Finally, the compatibility demand (|15p is reformulated in section 12.41 The resulting crite- 
rion can only be a necessary one since we deal with the problem on the level of instruments 
without referring to any experimentally realizable measuring process (cf. [31 1). 

2.1 Measurement outcome 

Given some measurement signal, say [Z]q, we must be capable of determining the corre- 
sponding trajectory of post-measurement states up to time t by assumption. So we assume 
the existence of a map 3>t from the set of all measurement records {[-Z]q} to the set of all 
normalized linear NMQSD solutions: r (->• |Vv([z*]o)) with r G [0, t] and l^oC^o)) = 
This map is surjective because of demand (ii). If two measurement signals [Zi]q, [Z^q are 
associated with the same trajectory of post-measurement states, the measurement signal 
can be post-processed such that only [Zi]g is displayed as the measurement outcome in 
both cases. Consequently, <3?i is a bijection w.l.o.g.. 

Any process realization [z*]q can be associated with a normalized solution of the linear 
NMQSD. Conversely, a given normalized linear NMQSD solution can be differentiated with 
respect to time. Due to the at most linear dependence of the ansatz operator on the noise, 
the corresponding process realization [z*]q can iteratively be determined by comparing 
that time-derivative with the known equation of motion. Hence, there is a bijection, say 
St, from the set of all normalized linear NMQSD solutions within [0,t] into the set of all 

2 The proof presented in this paper is more general: It also holds for ansatz operators which feature 
an additional linear functional of the noise. In the context of energy transport in molecules, a functional 
Taylor expansion in the noise has recently been applied in order to approximate the ansatz operator |16| . 
The line of argument presented here is therefore valid for systems for which this Taylor series terminates 
after the second term. 
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process realizations {[z*]o}- As a result, £t o <£ t provides a one to one correspondence 
between the measurement signals [Z] and the process realizations [z*] . For this reason, 
[z*]q can be assumed as the measurement signal w.l.o.g.. 

2.2 Stochastic evolution operators 

Let Gtdz*]^) be the so-called stochastic propagator which maps the initial system state 
|V>o) to the linear NMQSD solution at time t, i.e. \M[ z *]o)) = Gt([z*] )\ipo) . It is easy to 
realize that the ansatz operator 0(t, s, [z*] ) exists if and only if Gt([z*]o) is invertible: If 
Gt([ z *]o) is non-singular, the ansatz 

0(t,s,[z*}l) = (JL&([**]&)) ~* (21) 

certainly satisfies ©. Conversely, if 0(t, s,[z*] ) exists, one can integrate the evolution 
equation for G^Qz*]!)), i.e. effectively (J7|). So (-^([z*]^) turns out to be a time-ordered 
exponential, which is indeed invertible. This observation allows us to define a two times 
propagator for any < s < t 

^([z%):=dt([zr]l){d a ([^))-\ (22) 

The two times propagator obviously obeys A*([z*]q)|t/ ; s([2 : *]o)) = I^M[ z *]o)) an d satisfies 
the cocycle condition 

At([z%) = Al([z%) A r s ([z*Yo) (23) 

for any < s < r < t. With a glance at postulate (|18|) . we can rewrite demand (ii) 
equivalently as V^Qz*]^) oc ^4q([z*]q). Without loss of generality, the absolute square 
of the proportionality factor can be absorbed in the scalar measure Hq(-)- So up to an 
irrelevant phase factor, we may identify Vo([ z *]o) = j4q([z*]q). 

2.3 Scalar measures 

Recalling that Vq(-) denotes the path probability measure of the Gaussian process [z*]q, 
we explicate the minimality demand as follows: Any event Fq 6 Ft with Vq(Fq) = is 
verified in a measurement up to time t with the probability Pq(Fq\p) = 0, independently 
of the pre-measurement state p. So demand (iii) translates to an absolute continuity of 
the operator- valued measure (l" *)t (•)[!] with respect to Vq{-)'- 

VF^eFf. v&F*) = => (Y$)\f*)[1] = 0. (24) 

The finite, positive scalar measure p can be decomposed into a singular part /4 and an 
absolutely continuous part /i* with respect to v\ by means of Lebesgue's decomposition 
theorem. The singular part /i* is concentrated on a set Cq £ Ft with Vq{Cq) = and 
/4(Oq \ Cq) = 0. Due to ([M]), one finds for the dual instrument 

(y *) t (^)[l]= / (A ([z%)y A ([z%) ^(d[z%) = 0. (25) 
Jc* 

The integrand of (|25p is a strictly positive definite operator as the two times propagator is 
invertible. Therefore, /4(Co) = must hold and, consequently, /4( - ) vanishes everywhere. 
According to the theorem of Radon-Nikodym, we may conclude that there is a fQ-a.e. 
unique positive functional fo([z*} ) with /Ug(d[z*]o) = fo([ z *]o) u o(^l z *}o) ■ 

How can these functional /o([^*]o) be interpreted? Whereas V ([z*} ) = A ([z*]q) 
ensures that the post-measurement states depend on the measurement signal in the same 
way as the normalized linear NMQSD solutions depend on the process realization, ,/o([z*]q) 
determines the statistics of the measurements due to (fT7|) . Consequently, the family of 
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functionals /o([^*]o) can be seen as a means to implement the non-demolition property 
([8]). In addition to this, /o([<£*]o) can be regarded as the Jacobian of a transformation 
from the actually measurement signal [Z]q to the process [z*]q. However, not every choice 
of these functionals results in a consistent family of instruments, which is the issue of the 
next section. 

2.4 A necessary criterion 

According to the previous sections, the normalized linear NMQSD solutions can only be 
generated as post-measurement states if the instruments representing the time-continuous 
measurement are of the following form 

Y*(Ft)[p]= ! Ai(\z%)p(Ai([z%)y m*%Mw?]t>)- ( 26 ) 

Each map ()26j) is an instrument if and only if it is normalized, i.e. 

(Yt)\nl)[t}= [ {Ai([z%)y Ai([z*]l) f t o([z*] t )u t (d[z*]l)=t, (27) 

which is a constraint on /o([^*]o)- Moreover, the compatibility (|15p has to be ensured. Let 
us fix < s < t and some Fq G J s . Omitting all arguments, we rewrite the compatibility 
condition as 

= E v s[l F sE ui [f (J§) t Ai]]. (28) 

Here Ia(') refers to the indicator function of the set A and E„t[-] denotes the conditional 
expectation value given some path [z*]q. As (|28p holds for any Fq £ JF s , one may conclude 

/ S (ig) f i§ = E vi [f (i*o) f i*] ^-a.e.. (29) 

Because of the cocycle property ([23]) and the two times propagator being invertible, the 
last identity may be written as 

l/o s (re= / (AKlz*]^ Al([z%) /o*([*t)^(d[*tK) ^o-a-e., (30) 

where [z*]q) is the conditional path probability measure of [z*]q given a path realization 
[z*]q. To sum up, ()26p represents a time-continuous measurement scheme for every family 
of positive functionals /o([ z *]o) w hich obey the families of functional integral equations 

PZBDD - 

Now we tighten demand (iii) and require that the probability (density) to measure [z*]q 
shall coincide with the probability weight of the corresponding normalized linear NMQSD 
solution in the unravelling ([8|). This requirement necessarily leads to the choice /o([z*]o) = 
1. Then the normalization (|27p is fulfilled by construction |14| . In the Markovian regime 
a(t — s) = n5(t — s), the resulting family of instruments ([26]) coincides with that one 
presented in |20| . In this limit, Barchielli and Gregoratti have proven in [20] that the 
compatibility is fulfilled for any Hamiltonian H and any coupling operator L. 



10 



3 Application of the criterion to concrete systems 



In the following, we check the compatibility demand for /q([^*]o) = 1 in the non-Markovian 
regime. Therefore, we study two simple models for which all calculations can be performed 
analytically without any approximations. 

3.1 Jaynes-Cummings model 

The original Jaynes-Cummings model describes the interaction of a two level atom with 
a single electromagnetic mode in the dipole interaction and rotating wave approximation 
|36j . Here we consider a generalized Jaynes-Cummings model where there is a coupling 
of the atom to arbitrarily many environmental modes. In particular, this model has been 
studied by means of NMQSDs [14]. For this problem, the system Hamiltonian reads 
H = u>/2a z , where u is the atomic transition frequency and o z the Pauli spin matrix of 
the z-direction. The coupling operator is given by the lowering operator L = <r_. Then the 
ansatz operator also turns out to be noise independent: 0(t,s, [z*]q) = /(i, s)<7_. Here 
f(t, s) is a function obeying /(s, s) = 1 and 

d t f(t,s) = (ico + F(t))f(t,s), (31) 
where F(t) := f ds a(t — s) f(t, s). The corresponding linear NMQSD equals 

dt\Mlz%)) = (-&** + A »- - Ht) *+*-) (32) 



Due to the noise-independence of 0(t, s,[z*]q), the master equation can be determined 
from ([32]) [37]: 

dtPred{t) = - i [ ^<5- 2 + 3m(F(£))(T + <T_, p ied (t)) (33) 
+ Wt(F(t))([&-p tei (t), &+] + [*_, p Ied (t)a+\). 

Consequently, the evolution equation for p Te d(t) is of the Lindblad form but with a time- 
dependent Lindbladian in general. Only for a(t — s) = n5{t — s), the Lindbladian becomes 
time-independent because one has then F(t) = k/2. For a non-Markovian open quantum 
system, the decay rate D\e(F(t)) can become negative, which does not, however, cause any 
problems in the NMQSD approach. In this context, see also [12J. 

One easily finds the two times propagator for the Jaynes-Cummings model, namely 

-t * t ( e -if (*- s )-/ s tdTi? ( T ) \ 

A S ([Z ]o) = ( 4-iHf.s-Ml rt , * -v,,t- f T At' F(t'\ 4-iH( t _ s ) ) > ( 34 ) 



e +if (-+*) Jt dr z * e -iwr-J7 dr ' F ( r ') e +i 



where (1, 0) T and (0, 1) T represent the eigenvectors of a z corresponding to the eigenvalues 
+1 and —1, respectively. Remarkably, the two times propagator depends only on the 
noise realization within [s,t]. This property has been identified by Barndorff-Nielsen and 
Loubenets as being one of two conditions which lead to an instrumental process with 
independent increments if they are simultaneously fulfilled - given that the measurement 
scheme is consistent, of course [34] . The integrand of the compatibility demand turns out 
to be of the form 

(ii«)) + i^To) = ( $$$y ) • (35) 

The functionals fo*(-) and j*(-) are stated in the appendix. According to the previous 
section, the following conditions must be satisfied if the normalized linear NMQSD so- 
lutions allow a measurement interpretation: E^t [/i*] = 1 and K u t [j l s ] = for any given 
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realization [z*]q. We have explicated both conditions by discretizing the occurring path 
integrals and taking the left boundary points of the tiny time intervals as intermediate 
points. The resulting high-dimensional Gaussian integrals can be evaluated analytically. 
Having performed the continuum limit again, the condition E„t [j*] = turns into the 
following integral equation for the correlation function 

dr a(r - u) exp ( iur - / dr' F* (r) ) = (36) 

J s 

for any < s < t and any u € [0, s). This integral equation is certainly satisfied by the 
correlation function of a Markovian open quantum system, i.e. a(t — s) = n5{t — s). For 
this special case, an easy calculation verifies E^tf/t*] = 1. So the normalized solutions 
of (|32p allow a measurement interpretation in the Markovian regime as expected (cf. 
[3"ll4"tf20] ) , If, however, we restrict ourselves to functions a(-) which are Lebesgue integrable, 
the fundamental theorem of calculus tells us that the only solution of ([36]) is a(r) = 
Lebesgue-a.e.. In this sense, the example shows that the normalized linear NMQSD 
solutions generically allow a measurement interpretation only in the Markovian regime. 



3.2 Dephasing interaction 

Let us consider one of the simplest open quantum systems in order to analyse whether 
there are certain limits in which a measurement interpretation is approximately possible. 
A qubit, described by H = uj/2 a z , shall be coupled to the bosonic environment in a purely 
dephasing way: L = r/ncr z oc H. Here k denotes some time scale and the dimensionless 
parameter r controls the interaction strength. It is easy to show that the ansatz operator 
is noise-independent for this setting: 0(t,s, [z*]q) = L [H]. One then obtains 

OtMVU = ( - £ + A r -?* ~ -2 f ds «(* - *)) hk([*1S)> (37) 

for the linear NMQSD, which can directly be integrated. Similarly to the previous section, 
the master equation can be stated and is of the Lindblad form with an in general time- 
dependent Lindbladian: 

2 

dtPred{t) = -i [-Z&z, Pred(t)] + —^{K(t)) ( [a z /3 red (t) , U x \ + [& z , p IeA (t ) U z ] ) , (38) 

where K(t) := f£ ds a(t — s). One finds again that the two times propagator depends only 
on the noise realization within [s,t]: 

2 ft 

A t s ([z*] t s ) = exp(-£(t-s)a z -^e(t,s) + - / drz*a z ), (39) 

Z K K J s 

where Q(t, s) := dr' J Q T dr" a(r' — r"). The discretized path integral of ([30j) with /q = 1 
can be evaluated analytically. The resulting effective compatibility demand reads 

V0 < s < t, VuG [0, s) : f dr a(r - u) = (40) 

J s 

in the continuum limit. So again a measurement interpretation is exclusively possible in 
the Markovian regime. 

It is nevertheless instructive to calculate the right hand side of the compatibility de- 
mand for the correlation function of an Ornstein-Uhlenbeck process, 

a(t-s) = fe-^s\, (41) 
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which tends to k 5(t — s) as T — > oo. If one divides (s, t] into N tiny intervals of the length 
At = (t — s)/N and defines z*- := z* + -^ v j € {0, N}, the conditional path probability 
measure can be expressed as 

N , 2 At N - 2 

= jzJli^m) «* (- £ E + M ) (42) 

i=i i=i 

With respect to this measure, the right hand side of the compatibility demand, say 
CnHs(s,t, [z*]q), turns out to be a path integral with a quadratic action. Such an expres- 
sion can be evaluated by employing standard techniques of path integration (e.g. [38]). 
which leads to 

C RHS (s,t, [z*] s o) = exp{ - ^[r\e- rt - e~ rs ) (43) 

2 

- 2r(a z Ke(z s ) - r)(l - e~ r ^) - j(l - e~ 2r ^)}}. 

As expected for a Markov process like the Ornstein-Uhlenbeck process, this expression 
depends exclusively on (the real part of) z* and not on the whole path [z*]q. The com- 
patibility is clearly fulfilled in the Markovian limit: 

lim C RHS (s,t,z* s ) = l. (44) 
r— s>oo 

Trivially, the compatibility holds in the limit of vanishing interaction r — > 0. The system 
dynamics then becomes unitary and independent of the measurement signal [z*]q, which 
means that the dynamics is induced by a trivial measurement. Moreover, the compatibility 
can be achieved approximately either in the weak coupling and short correlation time limit 
fir > 1, which is the physical Markov limit in fact, or in the short-time limit (t — s) <C 1/T 
and (t — s) -C k. In the latter scenario, the correlation function is approximately constant 
on the considered time scales, a(t — s) kT/2, which results in linear dependence of zt 
on z* almost surely, namely Zt ~ z* pa 0. By inspecting the two times propagator ([39]) . 
one can easily show that the system dynamics turns out to be effectively unitary: The 
considered times are much smaller than the time scale of the interaction. Therefore, the 
system may be regarded as being decoupled from its environment. Finally by inspecting 
(143 p . one realizes that there is no meaningful long-time limit such as (t — s) 3> 1/T which 
approximately leads to compatibility. 

For a fair discussion of these results, one has to pay respect to the involved time 
scales. Since the measurement is described by a family of instruments y, the continuum 
limit is anticipated. Consequently, the time span between two consecutive measurements, 
which constitute the time-continuous measurement scheme, is arbitrarily short and thus 
the shortest involved time scale - unless the open system is Markovian. In the latter case, 
the time span between two consecutive measurements is comparable with the time scale 
of the environmental dynamics. 

4 On the nonlinear NMQSD 

In this section, we show that the necessary criterion (|26l27|30p can also be used for inves- 
tigating a possible measurement interpretation of the nonlinear NMQSD in principle. 

Due to the linear dependence of the ansatz operator on the noise and the shift rela- 
tion (|10p . there is a bijection between the process realizations [z*]q and solutions of the 
nonlinear NMQSD up to time t. So [z*]q can be assumed as the measurement signal 
w.l.o.g.. Moreover, © indicates that |?/ ! t([5*]o)) depends on the realization [z*]q in the 
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same way as |V ; t([-2*]o)) depends on [z*]q. Therefore, we may conclude Vg (•) = -Ag(-) with 
a glance at (|18p . So the normalized linear NMQSD solutions can be turned into nonlinear 
NMQSD solutions by appropriately weighting the driving process realizations [z*]q. When 
reviewing the derivation of the nonlinear NMQSD [TJ], one realizes that the path proba- 
bility measure of the shifted process z% is absolutely continuous with respect to the path 
probability measure of z%, i.e. ^g(')- Consequently, the minimality demand (hi) implies 
/^o(d[-z*]o) = /o([ z *]o) zy o(d[- z *]o) a S a i n - As a result, any family of instruments which might 
generate the nonlinear NMQSD solutions as trajectories of post-measurement states must 
be of the form ([26 p . 

So if the nonlinear NMQSD can be interpreted in terms of a time-continuous measure- 
ment scheme, there must be a family of functionals /g([z*]g) such that (|27|) and (|30p hold 
and that the pure post-measurement state trajectories occur with the same probability 
(density) as they appear as nonlinear NMQSD solutions. 

5 Conclusions 

In this paper, we have considered NMQSDs whose ansatz operators depend at most lin- 
early on the stochastic process z\. For this class of open quantum systems, we have derived 
a necessary criterion for the measurement interpretation of non-Markovian quantum tra- 
jectories. By employing a representation theorem for instruments, any time-continuous 
measurement scheme which is conceivable in the framework of standard quantum me- 
chanics has been taken into account. In particular, this criterion covers both "direct" 
and "indirect" measurements. Here the attribute "direct" refers to general measurements 
which are implemented by some measurement apparatus in the sense of the dilation the- 
orem (cf. [32] ). In contrast to this, "indirect" measurements bear on settings where the 
given harmonic environment acts as a quantum probe (e.g. heterodyne detection). By 
tracing out the environmental degrees of freedom, any instrument on the bipartite Hilbert 
space of the "indirect" measurement setting reduces to an instrument acting on the system 
Hilbert space as long as the initial state is not entangled. Hence, the above formalism also 
applies in this case. 

A word of caution is in order here: The proof presented in this paper is based on the 
assumption that any measurement can be described by a completely positive instrument. 
Yet an entangled initial state of the total system could lead to indirect measurements upon 
the system which cannot be described by complete positive maps (cf. |39|). 

We have especially concentrated on the interpretation of the linear NMQSD: For a 
suitable choice of /q, the family of instruments (126 p generates the normalized solutions 
of the linear NMQSD. In the Markovian regime, the resulting measurement scheme coin- 
cides with the family of instruments presented in [2U], where the compatibility of these 
instruments is shown for any hnite dimensional system. 

By investigating examples, we have shown that in general the normalized solutions 
|^t([^*]o)} °f the linear NMQSD allow a measurement interpretation only if a(t — s) = 
K5{t — s). The compatibility of these measurements is violated in the non-Markovian 
regime. As our analysis is performed on the level of instruments without considering a 
concrete measuring processes, a descriptive explanation for this failure cannot be given. 
We may only conclude that in general causality is violated by any measurement scheme 
which might generate the normalized linear NMQSD solutions. At least, one can state 
certain limits in which the compatibility is approximately satisfied. 

Up to now, it remains an open question whether the nonlinear NMQSD solutions allow 
a measurement interpretation in the non-Markovian regime. We hope that (|26I27|3(J|) can 
explicitly be applied to some system in order to investigate this question systematically. 
Our investigation moreover shows a mathematical peculiarity of the unravelling ([8]) in 
comparison to (|lip in the non-Markovian regime: In section [3[ we have effectively proven 
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that ||V ; t([ z *]c))ll 2 is n °t in general a (.F^-martingale if a(t — s) ^ n8{t — s). This implies 
that the family of probability measures {/ ||V't([ z *]o)ll 2 ^o(d[ 2: *]o)>^ > 0} is not consistent. 
Consequently, one can numerically perform a Monte Carlo simulation of /5 re d(i) based on 
the normalized linear NMQSD solutions IV'tG- 2 *]^)) but there is no overall path probability 
measure for the driving process z^ in this case. The Girsanov transformation, however, 
shifts this inconsistency of the measures into the dynamics of the nonlinear NMQSD 
solutions. The effectively driving process z% is then Gaussian and one can show that 
Kolmogorovs extension theorem is applicable for any a(t — s). This difference in the 
unravellings does not occur for the Markovian quantum state diffusion as ||V't([- 2; *]o)ll 2 
turns out to be a (J r t)-martingale (cf. [20 1). 

It would be interesting to invert the question discussed in this paper. Namely, is there 
a family of functionals /q such that the functional integral equations (|27l30p are satisfied 
and the a priori state Yq (Qq) [|"0o XV'o |] coincides with /S re d(t)? 

Finally, it is important to notice that the proof of (|26l27l30p mainly bases on measure- 
theoretical arguments. Therefore, the criterion (|27l30p can be applied to a broader class 
of stochastic Schrodinger equations, which is defined by the following two demands: First, 
any of these stochastic Schrodinger equations is either linear or can be formulated in a 
linear version. This linearity implies that the corresponding stochastic propagator is in- 
vertible. Second, this linear equation must be linearly driven by some stochastic process. 
For any of these stochastic Schrodinger equations, (|26|27|30p provides a tool to investi- 
gate whether the randomness of its solutions possibly originates from a time-continuous 
measurement. 
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Appendix 

For the Jaynes-Cummings model, the operator-valued integrand of (|30p contains the fol- 
lowing two functionals 
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